On the role of the symmetry parameter /3 in the strongly localized regime 
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The generalization of the Dorokhov-Mello-Pereyra-Kumar equation for the description of trans- 
port in strongly disordered systems replaces the symmetry parameter j3 by a new parameter 7, 
which decreases to zero when the disorder strength increases. We show numerically that although 
the value of 7 strongly influences the statistical properties of transport parameters A and of the 
energy level statistics, the form of their distributions always depends on the symmetry parameter /3 
even in the limit of strong disorder. In particular, the probability distribution is p(A) ~ A 13 when 
A — > and p(A) ~ exp(— cA 2 ) in the limit A — ► 00. 

PACS numbers: 73.23.-b, 71.30.+h, 72.10.-d 



It has recently been shown by Muttalib et al**^ 
that the Dorokhov-Mello-Pereyra-Kumar Equation 3 -^ 
(DMPKE) for the description of electronic transport in 
disordered quasi-one-dimensional systems can be general- 
ized to comprise also the strongly disordered case. They 
found that the main difference between the diffusive and 
localized regime is reflected in the spatial distribution of 
the electrons inside the sample. The DMPKE was de- 
rived under the assumption that the electron density is 
homogeneous. This assumption is valid in the limit of 
weak disorder (diffusive regime) and leads to universal 
behavior of the electron transport, which is determined 
by only three parameters: the ratio L z fl of the system 
length L z to the mean free path £, the number of scatter- 
ing channels N, and the symmetry parameter j3 of the re- 
spective random matrix ensemble. The latter determines 
the statistical properties of the the model, for instance, 
the fluctuation of the conductance, var g ~ /3 _1 . 

The derivation of the DMPKE is based on the following 
representation of the transfer matrix^ which describes 
the scattering of electrons coming from the left (right) 
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where A is a ./V-dimcnsional real diagonal matrix and the 
structure of the matrices u is given by the physical sym- 
metry. The diagonal elements A a (a = 1,2,..., N) de- 
fine the conductance of the system via the Economou- 
Soukoulis formula^ 
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The probability distribution p(X) of parameters A a is 
given by the DMPKEM 
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with J = Ila<6 l^a — A(,|' 3 . The parameterization A a = 
(cosha; a — l)/2 introduces a new set of variables x a 



{x\ < x 2 < • • •), which follows the Wigner-Dyson statis- 
tics. The probability distribution p(A a ) of the normal- 
ized differences A a = (x a+ i — x a )/(x a +i — x a ) is well 
described by the Wigner distribution^ 



Pp(0 = Ap^ exp(-BpC), 



(4) 



where for /3 = 1, 2, 4: Ai = tt/2, B x = tt/4, A 2 = 32/tt 2 , 
B 2 = 4/tt, A 4 = 2 18 /3 6 tt 3 , and B 4 = 64/9tt, respectively. 

Several investigations showed^— that the same func- 
tion also describes the probability distribution p(s) of 
the energy level statistics in the diffusive regime. Here, 
s is the difference of consecutive energy eigenvalues s = 
|ej+i — £i\/s divided by the mean level spacing s. 

In strongly disordered samples, the propagation of the 
electron is not diffusive. We cannot expect that all paths 
across the sample are equivalent. Mathematically, this 
leads to the re-formulation of the DMPKE into the more 
general form 



d(L,/£) J 



T ^ dX a 



X a (l + X a )K aa J 



dX a 



, (5) 



where parameters K a b depend on the statistical proper- 
ties of matrices u in Eq. ([T]). The explicit form of K a b is 
determined by the model symmetry^ - — The Jacobian 
J now have a form 
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Although Eq. (|5|) was derived only for orthogonal systems 
(j3 = 1), it can be shown to be valid also for (3 = 2 and 
/? = 4. The conductance is still given by Eq. ([2]), it 
becomes implicitly a function of the spatial distribution 
of the electrons. 

The main difference between the DMPKE and its gen- 
eralized version lies in the presence of the parameters 7 Q f, 
in the Jacobian. Later work showed^ that it is possible 
to approximate all 7 Q (, by a single parameter 7. Simi- 
larly, the parameters K aa are substituted by a constant, 
which is of order of unity in the limit of strong disor- 
der. It was argued^ that 7 — > /3 in the diffusive regime 
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but 7 —J- when the disorder increases. This assumption 
was confirmed, at least for the orthogonal symmetry, by 
numerical work.— 

If 7 really decreased to zero, one would expect that 
the probability distribution of A a and that of the level 
statistics should converge to the Poisson distribution 



Pp(C) 



(7) 



Such changes of the distributions due to the increase of 
the disorder are really observed both in the parameters 
Acr^ and in the level statistics. In the latter case it has 
been used for the estimation of critical parameters of the 
metal-insulator transitions ! 10 ' 15 ~ — 

The role of 7 and its relation to the symmetry parame- 
ter P still requires a more detailed discussion. Therefore, 
we investigate in this paper the shape of numerically ob- 
tained distributions p{C) in the limit of strong disorder. 
We show that, although 7 indeed decreases to zero, both 
the distribution of A = (x 2 — x\)/(x2 — xi) and the level 
statistics p(s) do depend on the symmetry parameter f}. 
In particular, the small £ behavior of these distributions 
is always given by 



p(0 ~ C p 



(8) 



independently on the strength of the disorder. However, 
this power-law behavior is observed only in a very narrow 
range close to zero. 

The systems to be investigated are defined on a 2D 
square lattice with lattice constant ai and described by a 
tight-binding lattice Hamiltonian with nearest-neighbor 
hopping terms 



y ] trr' c rtr c r'cr'- (9) 



Here, a = ±1/2 is the electron spin, r are the sites of the 
2D lattice of size L 2 , e r are the appertaining random on- 
site energies distributed according to the box probability 
Pbox^t) = (l/W)@{W/2 - |e r |), and W measures the 
strength of the disorder. Periodic boundary conditions 
are applied in both directions. For the symplectic Ando 



model, t™, is a 2 x 2 matrix with 
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and V? 



V? 



V 2 = 1. For the orthogonal model, 
t rr i = t°°, 8 aa i . Energies and lengths are measured in 
units of V and lattice constant ai , respectively. For E = 
and V\ = 0.5 the symplectic model exhibits a metal- 
insulator transition at a critical disorder W c ps 5.84.— 

The limiting behavior of the distribution p(() is bet- 
ter visible when the distribution of the logarithm of £, 
p(ln£), is studied instead. From the equation 

p(C)d( = p(Od£, £ = lnC (11) 
we obtain that the the relation p(() ~ £^ corresponds to 

mo ~ (fi + % e^-oo. (12) 
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FIG. 1: (Color online) The distribution p(A) for the 2D Ando 
model in the metallic, critical (W = 5.84), and localized 
regime. The linear size of the system is L = 14. The data 
are compared with the Wigner surmise -Ri(A) and with the 
Poisson distribution Pp(A). While p(A) is very similar to the 
Wigner surmise in the metallic regime (W = 2), it resembles 
the Poisson distribution when the disorder increases. How- 
ever, for any disorder strength, p(A) decreases to zero when 
A -> 0. 



Similarly, the large-£ tail of the distribution can be ana- 
lyzed from the function 



ln[-lnp(£)] = a£ 



(13) 



with a = 2 and 1 for the Wigner and Poisson distribution, 
respectively. 

In the limit of strong disorder, ^ > 1, the typical 
conductance is given by the smallest parameter x\ as 
g « e~ Xl . The parameter x\ determines the localiza- 
tion length A as x\ = 2L/X (L ^> A). Thanks to this 
relation, the transport properties of strongly disordered 
system can be understood from the numerical analysis of 
relatively small samples, provided that L ^> A. 

We analyze statistical ensembles of iVstat ~ 10 s square 
samples^ (typical size is L = 14) and collect the statisti- 
cal distribution of the normalized difference. The results 
are displayed in Figs. [T]|3] for the 2D Ando model and 
for the 2D orthogonal model. Fig. Q] exhibits the distri- 
bution p(A) for various strengths of the disorder. The 
data show that for small disorder (W = 2) the distribu- 
tion is very similar to the Wigner surmise. Although the 
form of the distribution changes when disorder increases, 
the decrease p(A) —> is noticeable even for W = 50. 
The small-A behavior of the distribution is better visi- 
ble in Fig. [5] which plots the distribution p(lnA). Our 
numerical data for any disorder show that the distribu- 
tion p(\n A) becomes parallel to the Wigner surmise for 
very small A. This proves that lnp(ln A) ~ (j3 + 1) In A 
and, consequently, p(A) ~ A' 3 (Eq. (fT2"j)). However, the 
power-like behavior p(A) ~ A 13 is observed only for a 
very small part of the statistical ensemble. For instance, 
the linear behavior p(\a A) ~ (f3 + 1) In A is observable 
only for In A < lnA m = -4 (W = 10), -7 (W = 20) 
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FIG. 2: (Color online) The distribution p(lnA) for the 2D 
Ando model. The linear size of the samples is L — 14. The 
solid curves represent the Wigner P4 and the Poisson Pp 
distributions. With increasing disorder, the distribution of 
p(lnA) changes. It is similar to Poissonian in the bulk, but 
for very small A we see the linear behavior lnp(ln A) oc In A. 
The straight solid lines are fits for W = 5.84, 10, 14, and 20 
with slopes 4.789, 3.78, 3.0, and 3.3, respectively. 



(Fig. [2]). The probability p to find a system with such 
a small value of A decreases rapidly when the disorder 
increases: p ~ lfT 3 (W = 10) but p ~ 1CT 6 for W = 20. 

According to single parameter scaling theory^ the 
same change of the distribution function is expected 
when the size L of the system increases while disorder 
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FIG. 4: (Color online) The same as in Fig. [T] but on a loga- 
rithmic scale for various strength of the disorder and L — 14. 
Solid lines are Wigner and Poisson distributions. For the lo- 
calized regime (W > 10), the mean value of the logarithm of 
the conductance is given in the legend. 



is fixed. Owing to the necessity to analyze huge statis- 
tical ensembles, we did not study the size dependence of 
p(A) for the Ando model. However, we checked the in- 
dependence for 2D orthogonal systems (Fig. [3]) . Again, 
the distribution p(lnA) follows Eq. (fl~2"|) with /3 = 1 pro- 
vided that A is sufficiently small. Also, our data support 
the scaling idea: two distributions are similar if they cor- 
respond to systems with the same value of (In g) . 

To estimate the large-A form of the distribution, we 
plot in Fig. |4]p(A) for the Ando model and compare it 
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FIG. 3: (Color online) The distribution p(ln A) for a 2D or- 
thogonal model with anisotropic hopping (t xx i = 0.9t zz /). 
Two ensembles with different values of L and W but with 
the same value of (lng) possess the same distribution p(A). 
The Wigner surmise Pi and the Poisson distribution Pp 
are also plotted. Fits shown by straight lines confirm that 
lnp(lnA) oc In A with slopes given in the legend. There- 
fore, p(A) ~ A when A — ► 0. The number of samples is 
iVstat = 1.6 x 10 9 for L = 20 (W = 12), and about 10 8 else. 



FIG. 5: (Color online) The energy level statistics p(lns) of a 
strongly disordered Ando model. In the limit of s — > 0, a clear 
~ s l3+1 behavior is observed for W = 10 from the comparison 
with the Wigner surmise P4 for symplectic symmetry. With 
increasing W the range where this agreement holds shifts to 
smaller a. The Poisson fit (upper solid line) is valid only 
for larger s. The system size is L — 20 and the disorder 
strengths shown are W = 10 (+), 13 (x), 15 (*), and 25 (□), 
respectively. 
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FIG. 6: (Color online) Plot of parameters xi and X2 for a 
statistical ensemble of 10 4 samples. The size is L = 200 and 
disorder W — 8. The inset shows the value of the difference 
X2 — Xi as a function of xi. Solid lines mark mean values, 
(xi) — 14.09 and {x2 — xi) = 6.08. The data confirm the level 
repulsion since X2 and xi never coincide. Note that small 
values of A are observed only in samples with relatively large 
value of x\ > (xi). Such samples only marginally influence 
the mean conductance. 

with the Wigner distribution P4 and the Poisson distribu- 
tion. Again, our data confirm that p(A) is never identical 
with the Poisson distribution. For large values of A the 
distribution is p(A) ~ cxp(— cA b ) with exponent b « 2, 
at least in the limit of A ^> 1. 

The numerical investigation of the energy level statis- 
tics generated a similar result. For large disorder, W > 
W c ~ 5.84, the large- s part of the level statistics p(ln s) 
is well described by the Poisson distribution as shown in 
Fig. O In the opposite limit s 0, a behavior close to 
the Wigner surmise P4 is observed in which the range of 
the agreement is continuously diminishing with increas- 
ing disorder. For very large W > 20, however, only the 
downturn can still be noticed. The eigenvalues have been 
calculated within an energy interval —0.5 < £j < 0.5 by 
direct diagonalization of the respective (2L) 2 x (2L) 2 ma- 
trices with up to 3 x 10 5 realizations. Additional calcu- 



lations for larger system sizes corroborated the results 
shown here for L = 20. 

Our numerical data indicate that the generalized 
DMPKE fails to describe correctly the small-A behavior 
of p(A). Contrary to the expected behavior p(A) cx A 7 , 
we observe for any disorder that p(A) <x A^. This re- 
striction influences the statistical properties of the con- 
ductance only weakly since samples with small A rep- 
resent only a very small part of the total statistical en- 
semble of the AT st at samples. Also, as shown in Fig. H3 
samples with small A possess a large value of the small- 
est exponent x\ > (xi), and consequently have also a 
small conductance. Therefore, their contribution to the 
conductance statistics is negligible. We expect that the 
GDMPKE gives a correct value of the mean In g but its 
description of the small-g-tail of the distribution p(hi g) 
can eventually differ from numerical (experimental) data. 

The aim of this paper was to verify that the physi- 
cal symmetry of the model governs the small £ behavior 
of the distribution p((), where C denotes the normalized 
differences A = a; 2 — X\j(xi — x\) or, for the energy 
level statistics, the normalized difference of consecutive 
eigenvalues s = — £i\/s. This was done by a nu- 

merical study for the metallic and also in the strongly lo- 
calized regime. Our results confirm that the generalized 
DMPK equation is not in contradiction with conclusions 
provided by random matrix theory. 

We also showed that the distribution p(A) never cor- 
responds entirely to the Poisson distribution, although 
7 — > 0. p(A) is not universal in the strongly disordered 
limit since 7 depends on the disorder. For small A wc 
found a distribution p(A) ~ A* 3 and in the limit of large 
A it behaves as exp(— cA 2 ). This observation is also con- 
sistent with the generalized DMPK equation. The Pois- 
son distribution indicates that the two parameters, x\ 
and X2-, are statistically independent. On the contrary, 
as discussed in previous work ; 12 ' 21 the statistical correla- 
tions survive for any disorder strength and are responsi- 
ble for the non-Gaussian distribution of the logarithm of 
the conductance. 
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